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Abstract: 

A two-dimensional array of quantum dots in a magnetic field is considered. The electrons 
in the quantum dots are described as unitary random matrix ensembles. The strength of 
the magnetic field is such that there is half a flux quantum per plaquette. This model 
exhibits the Integer Quantum Hall Effect. For iV electronic states per quantum dot the 
limit N — > oo can be solved by a saddle point integration of a supersymmetric field theory. 
The effect of level statistics on the density of states and the Hall conductivity is compared 
with the effect of temperature fluctuations. 

PACS Nos.: 71.55Jv, 73.20Dx, 73.40Hm 

We consider a two-dimensional array of quantum dots in a homogeneous magnetic field 
perpendicular to the array. A quantum dot in an array is a complex finite system of 
electrons, subject to strong Coulomb interaction and a confining potential. Even if the 
number of electrons is small there is a large number of electronic states in a given energy 
interval. Therefore, we are forced to use a statistical description of the quantum dot. A 
typical feature of such a complicated non-integrable system is level repulsion. The latter, 
also found in other complex many-particle systems like atomic nuclei [1], atoms [2] or 
metallic particles [3] , can be conveniently described by random matrix ensembles [4] . Since 
the magnetic field breaks the time-reversal invariance in the dot, an appropriate model is 
the Gaussian unitary ensemble (GUE). Electrons can travel in the array of quantum dots 
due to tunneling between neighboring dots. On the square-array, which will be considered 
in this article, the tunneling rates are t and t' for nearest and next nearest neighbors, 
respectively (cf. Fig.l). The coupling between the individual quantum dots due to these 
tunneling processes is weak. This allows us to assume that the statistical occupation of the 
electronic states in each dot is uncorrelated between different dots. Thus the quantum dots 
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can be represented by independent random matrix ensembles. Moreover, we also assume 
for simplicity that the tunneling processes are independent, i.e., the tunneling electrons do 
not interact with each other. 

For very weak tunneling rates the array should behave like an insulator because of the 
fluctuations of the energy levels. One would expect for increasing tunneling rates that a 
metallic regime can be reached where the array becomes conducting. However, due to the 
statistical fluctuations of the energy levels in the dots the effect of Anderson localization 
must play a crucial role in the array. Anderson localization prevents a two-dimensional 
system to become metallic, at least if no or only a weak magnetic field is present [5] . On the 
other hand, in the two-dimensional electron gas in a homogeneous magnetic field quantum 
Hall transitions (QHT) have been observed which are accompanied by delocalized electronic 
states [6]. A QHT occurs if a gap opens in a band of electronic states. This phenomenon 
is known, for instance, from electrons which are subject to a homogeneous magnetic field 
and a periodic potential [7]. Depending on the magnetic field the electrons form several 
subbands where each subband contributes e 2 /h to the Hall conductivity [8,9,10]. As an 
approximation of the periodic potential one can use a tight binding model where the lattice 
constant is given by the period of the potential. In this article we will study the effect of 
the statistics of energy levels and the effect of thermal fluctuations on the QHT. 

There are two different approaches to the transport in quantum dots. One is based on 
the S-matrix, the other one on the Hamiltonian. The former is very useful for numerical 
simulations because it describes directly the reflection and transmission of the electrons 
through the quantum dots [11,12]. The latter, however, requires the application of linear 
response theory to get a conductivity via Kubo's formula. In this article the Hamiltonian 
representation will be used. The effective Hamiltonian of an array of quantum dots reads 
as a quadratic form ^2H r ' r , c^c^, in the fermion creation and annihilation operators c\ 
c with the matrix elements 

H^y = Hy a 5 r y + H' rr ,5 a ^ a > +V r 5 r y5 a ^ a > , (1) 

where a, a' = 1,...,N label the iV electronic states in the quantum dots and r and r' 
label positions of the quantum dots in the two-dimensional array. In general, tunneling 
between all N states should be allowed with some probability, depending exponentially 
on the energy of the states a and a'. To include this would require a detailed knowledge 
of these states. Therefore, we choose as a simplifying approximation the assumption that 
there is tunneling only between states with the same a at nearest or next nearest neighbor 
dots with fixed tunneling rates. The distance between neighboring dots is measured in 
units of (4> /2B) 1 ^ 2 . Typical distances are a = 100...500nm [13]. The magnetic field for 
the creation of one flux quantum per plaquette is B = (po/a 2 0.016...0.4T. This regime 
is accessible in natural crystals (a ps 0.5nm) only with astronomical magnetic fields. 
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The electron can occupy statistically states inside the quantum dot which are repre- 
sented by the matrix elements H^ ,a . H is the N x N Hermitean Hamiltonian (H^ = H) 
of a quantum dot with N 2 statistically independent matrix elements. They are Gaussian 
distributed with zero mean and {H^ a H^, a ) = (g/N)5 aa 5 a a 8 r y. g is the strength 
of the level fluctuations. This depends on the strength of the interaction beween the elec- 
trons inside the dot. Therefore, g increases with the number of electrons per dot and with 
increasing confinement. 

The tunneling is represented by the Hamiltonian H' . This reads in Landau gauge 
(with r = (x,y)) for flux per plaquette 

H' r r , = te 2 ^^°5 r ,, r+ex + t8 r ,, r+ey ± it'e 2i ^^5 r ^ r+ex±ey + h.c. (2) 

For the special case of half a flux quantum per plaquette ((f) = <j>o/2) the phase factor in 
(2) is real and changes only sign between nearest neighbors in y-direction. Finally, the 
potential term V r represents an additional external (e.g. electric) field. Here we regard a 
staggered chemical potential V r = (— l) x+y n which opens a gap 2fi in the spectrum of the 
electrons, as will be explained below [14] . It is probably difficult to implement a staggered 
potential in a real sample of quantum dots. However, the parameter \i plays only the role 
of a gap which could be created also by other means. 

We choose for the tunneling rate t = 1. Therefore, [i, t' are measured in units of t and 
g is measured in units of t 2 . If we identify fermions with the four corners of the unit cell 
(Fig.l) the tunneling matrix H' can be diagonalized by a Fourier transformation. This gives 
a 4 x 4 matrix in Fourier space. H, the Hamiltonian of a dot, is a diagonal matrix with re- 
spect to the four corners in the sublattice representation H = (H" ,a , H^'" , H^ ,a , H% ' a ). 
A similar model with correlated randomness H^ ,OL = H^' a = —H^'" = —H^' a was con- 
sidered in Ref. [15]. 

We begin the discussion of the model with the analysis of an array where the inter- 
action of the electrons inside the quantum dots are neglected. It can be understood as a 
tight-binding model for non-interacting electrons in a metall with some electronic bands in 
a magnetic field [14,16]. The Fourier components of H' can be expanded around the four 
nodes k = (±7r, ±7r) for k = (±7r, ±7r) + ap with small p vectors. After a global orthogonal 
transformation the Hamiltonian reads 



H"{p) = 2 



/ ii — t' ipx-Py -2t'(p x +p y ) \_ 

-iPx-Py -fJ> + t' -2t'(p x -p y ) 

-2t'(p x +p y ) fM + t' Py+iPx 
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The last equation combines the 4 x 4-structure to a 2 x 2-structure with 2x2 block matrices 
H'/,. Neglecting terms 0(p 2 ) the Green's function (H + iuj)~ 1 decays into a diagonal block 
structure 

G ^)-{ (H^l N + h l+ tu)^) (4) 

with the diagonal matrix h\ = (Hi + Hj >1 U\2 + H4). Thus the diagonal elements are 
statistically independent. Ijv is the N x A^-unit matrix. It is interesting to notice that the 
matrices H'-^ = mj<j z + z'V ' x u x =F 2V 7 y cr y represent two independent two-dimensional Dirac 
Hamiltonians with masses m\ = ji — t' and m<i = \i + t' , respectively. 

The current density in a Dirac model can be calculated from the response to an 
external vector potential q y [17]. The introduction of this vector potential is equivalent to a 
change of the boundary conditions, a concept extensively used in numerical investigations 
of Anderson localization [18]. The response to the vector potential leads to the Hall 
conductivity a xy in terms of Green's functions. We obtain for q y ~ the expression 
[14,15] 

Cx y ~ i / Tr[a x G rr i(E — iu)G r ' r "(E — iu)a y G r " r (E — zu;)]— . (5) 

r,r',r" 

Here ^ s the sum normalized with the number of quantum dots in the array and the 
number of energy levels N. If there is only one electron per dot the energy spectrum has 
discrete levels which are well-separated. For instance, with a harmonic oscillator potential 
for the dot we have E n = fko p (n + 1/2). The separation of the energy levels in the single 
electron case allows us to neglect all levels with n > 0. Consequently, there is no statistics 
of energy levels and we can write h\ = 0. For the Hall conductivity in units of e 2 /h we 
find 

a xy = (l/2)[sign(mi)e(|mi| - \E\) + sign(m 2 )e(|m 2 | - \E\)}, (6) 

where G is the Heaviside step function. This result reflects correctly the qualitative be- 
havior of the Hall conductivity at the QHT: The Hall conductivity of the original lattice 
fermion problem is the sum of the Hall conductivities from the light Dirac mass (mi) and 
the heavy Dirac mass (777-2), such that the total a xy has a jump from to 1 if the light 
mass changes the sign (i.e., exchange of particles and holes in the Dirac model). Thus the 
Dirac fermions, together with the Hall conductivity of Eq. (5), represent a simple picture 
for a Hall transition. Special cases are \x = which gives a xy = and the (unrealistic) case 
t' = with a xy = (sign(//)/2-7r)6(|mi| — \E\). The sharp step-like QHT is only possible 
in an ideal systems of non-interacting lattice electrons at zero temperature. In order to 
compare with real systems we have to include the statistical fluctuations of the energy 
levels as well as thermal fluctuations. The latter are taken into account by replacing the 
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integral over u> in (5) by a summation over discrete Matsubara frequencies u n = (2n + l)7rT 
(n = 0, ±1, ±2, ...). This leads to a thermal broadening of the step-like behavior of a xy . 

The effect of the level fluctuations is evaluated by averaging a xy over the random 
matrix elements of h\. In order to derive a simple expression for the limit of infinitely 
many energy levels per dot (N — > oo) it is convenient to write the product of Green's 
functions G = (Hq + hi + zero) -1 (Hq is either or K'^) in the expression of the Hall 
conductivity formally as a functional integral of a supersymmetric model [15,19,20] 

G$G%,,(P* T = (KK'X^r'^J»K)s - (tfx^^x^x^s (7) 

with (...)s = / ... exp(— Si) Y[ r d& r d& r and with the supersymmetric action (sum conven- 
tion for a) 

S! = -tS z J2 *?,Mj(tf0 + **0)rj;r'j'^ (8) 

where s z = sign(Im^) and the field ■ = (\&" - •). The first component is Grassmann 
and the second complex, fj, = 1,2 labels the complex and the Grassmann components, and 
j = 1, 2 labels the two components of the Dirac model. This choice guarantees a normalized 
functional. Consequently, the averaging with respect to the Gaussian distributed matrix 
elements of h\ can be performed in the functional integral as (exp(— S\))h x = exp(— S2) 
with the effective action S 2 . The latter is obtained from Si by replacing the second 
term with (g/N) ^ r j(®r,n,j®?,n,j) 2 - Thus we have derived an effective field theory for 
$ which serves as a generating functional for the average product of Green's functions. 
It is important to notice that not only hi creates the interaction in S2 but also other 
types of random terms in Si. For instance, the interaction can also be created by a term 
(N/g)(Q r ^j) 2 — 2is z Q r ^ j j$", f i j^r,^ ,j as ^ ne second term in Si, followed by an integration 
over the matrix field Q. This field, in contrast to the random matrix hi, does not depend 
on the index a of the electronic states inside the quantum dot. This means that the 
distribution hi can be transformed into another distribution with a new 'random variable' 
Q (which does not have a probability measure but some generalized distribution including 
Grassmann variables). In other words, we can write, after integrating out the field $, 

([(Ho + hi+ zao)- 1 ]^-)^ = ([(H + 2Q + zao)' 1 }^ ...)q. (9) 

The distribution which belongs to (...)g was investigated in detail in [20]. Here we need 
only the result for leading order in iV: (...)q = J* ... exp(— NS(Q, P)) Y[ r dP r dQ r with 
diagonal matrix fields Q r , P r and 

S(Q,P) = - n ^ Tr ^ +Tr ( P r)] 
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+ log det(H + 2Q + za ) - log det(H - 2iP + za ) (10) 

The number of levels N appears in front of the action. Thus the effect of the statistics of 
the energy levels can be evaluated for N — > oo in saddle point (SP) approximation. The 
SP equation reads 

^-[-Tr(Q 2 r ) + logdet(H + 2Q + za )} =0. (11) 

A second SP equation appears from the variation of P by replacing Q — > — iP. As an 
ansatz we take a uniform SP solution Q = —iPq = (1/2) [17700 + M s cr 3 ]. Then (11) 
leads to the conditions 77 = (77 + ui — iE)gI, M s = —migl/(l + gl) with the integral 
1 = f[( mi + M S ) 2 + (r] + uj-iE) 2 + k 2 }- 1 d 2 k/2ir 2 . This result means that disorder shifts 
the frequency ui — > uj + rj and the Dirac mass mi — > M' = mi + M s , where 77(777-1,0;) 
and M s (mi, a;) are solutions of the SP equation. For instance, with w = we have 
rj 2 = (1/4) (M c 2 - m{)Q{M 2 - ml) where M c = 2e _7r / s . The sign of 77 is fixed by the 
condition that 77 must be analytic in uj. This implies sign(77) = sign(o;). The average 
density of states (DOS) is proportional to 77 in the — > oo-limit. Thus we have a narrow 
DOS for the array of quantum dots of width 2M C in contrast to the isolated dot which 
has a semicircular density of width 2^fg. The DOS vanishes for E = in the absence 
of level fluctuations. The creation of a non-zero DOS due to level fluctuations is a non- 
perturbative effect. 

At T = and E = the Hall conductivity per fermion level reads in the limit N — > 00 
and with the approximation that M' and 77 do not depend on u 

a xy Ri 1/2 + sign(mi ) [l/2 - (l/7r)arctan(^/M c 2 /m? - 1)0(M C 2 - m?)l . (12) 

The Hall conductivities are plotted in Fig. 2 for T = 0.1 with and without level fluctuations. 
It is remarkable that the Hall conductivity is enhanced by the level fluctuations for a xy < 
1/2 whereas it is suppressed for a xy > 1/2. The effect of these fluctuations is strictly 
constrained to the interval 2M C . 

Conclusions In a square-array of quantum dots with N electronic states per dot we 
have investigated the DOS and the Hall conductivity. Both quantities are significantly 
affected by the statistical fluctuations of the energy levels. In particular, the Hall conduc- 
tivity, which is step-like at the QHT in the absence of fluctuations, has a more complicated 
behavior in the presence of level fluctuations. Thermal fluctuations have a different effect 
on the Hall conductivity; they lead to a simple broadening of the step-like behavior. 

Only the average quantities have been considered. However, it is possible within the 
same method described in this article to study also higher moments of these quantities. 
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Figure Captions 



Fig.l: Schematic picture of an array of quantum dots with nearest (t) and next nearest 
neighbor (t') tunneling. The square denotes the unit cell of the translational invariant 
array with magnetic flux $ = <&o/2. 

Fig. 2: Hall conductivity a xy in units of e 2 /h as a function of the effective chemical potential 
m = ji — t' at temperature T = 0.1. The circles are without level fluctuations and the full 
curve is with level fluctuations with variance g = 1.36. 
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Fig.l: 
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